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Abstract 

Hule and Rio Cuarto are maar lakes located 1 1 and 18 km N of Poas volcano along a 27 km long fracture zone, in the Central 
Volcanic Range of Costa Rica. Both lakes are characterized by a stable thermic and chemical stratification and recently they 
were affected by fish killing events likely related to the uprising of deep anoxic waters to the surface caused by rollover 
phenomena. The vertical profiles of temperature, pH, redox potential, chemical and isotopic compositions of water and 
dissolved gases, as well as prokaryotic diversity estimated by DNA fingerprinting and massive 16S rRNA pyrosequencing 
along the water column of the two lakes, have highlighted that different bio-geochemical processes occur in these 
meromictic lakes. Although the two lakes host different bacterial and archaeal phylogenetic groups, water and gas 
chemistry in both lakes is controlled by the same prokaryotic functions, especially regarding the CO2-CH4 cycle. Addition of 
hydrothermal CO2 through the bottom of the lakes plays a fundamental priming role in developing a stable water 
stratification and fuelling anoxic bacterial and archaeal populations. Methanogens and methane oxidizers as well as 
autotrophic and heterotrophic aerobic bacteria responsible of organic carbon recycling resulted to be stratified with depth 
and strictly related to the chemical-physical conditions and availability of free oxygen, affecting both the CO2 and CH4 
chemical concentrations and their isotopic compositions along the water column. Hule and Ri'o Cuarto lakes were 
demonstrated to contain a CO2 (CH4, N2)-rich gas reservoir mainly controlled by the interactions occurring between 
geosphere and biosphere. Thus, we introduced the term of bio-activity volcanic lakes to distinguish these lakes, which have 
analogues worldwide (e.g. Kivu: D.R.C.-Rwanda; Albano, Monticchio and Averno: Italy; Pavin: France) from volcanic lakes only 
characterized by geogenic CO2 reservoir such as Nyos and IVlonoun (Cameroon). 



Citation: Cabassi J, Tassi F, Mapelli F, Borin S, Calabrese S, et al. (2014) Geosphere-Biosphere Interactions in 6/o-,Acf/wfy Volcanic Lakes: Evidences from Hule and 
Rio Cuarto (Costa Rica). PLoS ONE 9(7): el 02456. doi:l 0.1 371/journal.pone.Ol 02456 

Editor: Dwayne Ellas, Oak Ridge National Laboratory, United States of America 

Received April 2, 2014; Accepted June 19, 2014; Published July 24, 2014 

Copyright: © 2014 Cabassi et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper. 

Funding: This work benefitted by the financial support of the 7th Workshop of the Commission on Volcanic Lakes (lAVCEl), INGV and the Laboratory of Fluid and 
Rock Geochemistry of University of Florence (Resp. Franco Tassi). Francesca Mapelli was supported by University of Milan, DeFENS, European Social Found (FSE) 
and Regione Lombardia (grant "Dote Ricerca"). No funders were involved since this study was supported by using funds belonging to laboratories of different 
institutions, which participated to this research. Nevertheless, FSE and Regione Lombardia had no role in study design, data collection and analysis, decision to 
publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* Email: jacopo.cabassi@gmail.com 



Introduction 

Volcanic lakes are peculiar natural systems on Earth, although 
they are a common feature of volcanic systems characterized by 
recent activity, being present in 476 volcanic structures worldwide 
(VHub, CVL Group page; [1]). A volcanic lake simultaneously 
acts as both a calorimeter and a condenser for acidic volatiles from 
magmatic and hydrothermal degassing [2-6]. Thus, its existence 
and durability strictly depends on the balance between i) inputs of 
meteoric water and hydrothermal-magmatic fluids and ii) losses 



related to evaporation, permeation through sediinents and 
streaming [7]. Volcanic lakes were basically classified, as follows 
[1,4]: i) "high-activity" lakes affected by the addition of significant 
amounts of hot and hyperacidic hydrothermal-magmatic fluids; ii) 
"low-activity" lakes, characterized by C02-dominated fluid inputs 
at a relatively low rate from sub-lacustrine fluids discharges, 
favoring the establishment of a stable vertical stratification and 
possibly the accumulation of high amounts of dissolved gases in the 
deep water layers. At these conditions, a lake overturn triggered by 
either i) external events, such as earthquakes, landslides or extreme 
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weather conditions or ii) the progressive attainment of gas 
saturation conditions may cause the abrupt release of toxic gas 
clouds in the atmosphere. This phenomenon, also known as 
"limnic eruption", was firstiy documented at Monoun and Nyos 
lakes (Cameroon) in 1984 and 1986, respectively [8-15]. 
Accordingly, low activity lakes are commonly indicated as 
"Nyos-type" lakes. 

In Costa Rica, volcanic lakes are found in quiescent systems 
(Congo and Barva), as well as in volcanoes characterized by 
moderate hydrothermal activity (Irazu and Tenorio) and strong 
magmatic fluid emissions (Rineon de la Vieja and Poas) [16,17]. 
Hule and Rio Cuarto are low-activity, Nyos-type, maar lakes 
located at 11 and 18 km N of Poas volcano (Fig. 1), respectively, in 
relation of a 27 km long fracture zone passing through the Sabana 
Redonda cinder cones, the Poas summit craters (Botos, Active 
Crater and Von Frantzius) and the Congo stratocone [18]. In 
these two lakes, changes in the water color and fish death events 
were repeatedly reported, suggesting the occurrence of rollover 
episodes related to inputs of deep-originated gases [18]. To the 
best of our knowledge, no information is available on these lakes 
for chemical and isotopic compositions of dissolved gases deriving 
from geogenic sources and the structure of prokaryotic commu- 
nities. The latter are expected to play pivotal ecological functions, 
encompassing nutrient remineralization and carbon cycling, which 
is firmly linked to the fate of dissolved Ci gases, i.e. CH4 and COj. 

This paper presents the geochemical (water and dissolved gas 
chemistry) and microbiological results obtained from samples 
collected in 2010 during the 7* Workshop of the Commission on 
Volcanic Lakes (CVL; Costa Rica 10-21 March 2010), which is 
part of the International Association of Volcanology and 
Chemistry of the Earth's Interior (lAVCEI), by a group of 
geochemists, limnologists, biologists and volcanologists from 
different universities and scientific institutions. The aim of this 
multidisciplinary research was to unravel the bio-geochemical 
processes controlling the physical-chemical features of Hule and 
Rio Cuarto lakes along the vertical profiles, showing their 
implications for lake stratification and stability, and proposing 
evidences for a new classification system. 

Morphological and LImnological Outlines 

2.1 Morphological features 

Lake Hule (10°17'42"N, 84°12'37"W) lies within the 
2.3x1.8 km wide Hule basin, a volcanic depression also hosting 
Lake Congo to the north, which is separated from Lake Hule by a 
volcanic cone, and Lake Bosque Alegre (unofficial name) [18-20]. 
Lake Hule has a half-moon shape, a surface area of about 
5.5x10^ m^, an estimated water volume of 6.9x10^ m"*, and a 
maximum depth of ~23 m [17,18,21,22] (Fig. 2). The northern 
shoreline of the lake shows three tributaries, whereas an emissary 
(Rio Hule) is located to the NE [18,23,24]. 

Rio Cuarto maar (10°2r23"N, 84"13'OO"\\0 has a rim whose 
maximum elevation is ~52 m a.s.l. Lake Rio Cuarto shows steep 
sided walls and a flat bottom, a morphology typical of maar lakes. 
The lake has an E-W axis of 758 m, a mean width of 581 m, a 
surface of 3.3x10^ m^ and a water volume of 15x10® m^ [18,25] 
(Fig. 3). Rio Cuarto is the deepest (~67 m) natural lake in Costa 
Rica [19]. A small tributary is located on the eastern shore, 
whereas no emissaries were recognized [25]. 

The main morphological features of Hule and Rio Cuarto lakes 
can be summarized using the "depth-ratio" [26], which is a 
dimensionless parameter equal to the ratio between the average 
depth (the volume divided by the surface area of the lake) and the 
maximum depth of the lake. The obtained results are 0.55 and 



0.68, respectively, for Lake Hule and Lake Rio Cuarto, 
corresponding to an average depth of 12.6 and 45.5 m. According 
to Carpenter's heuristic classification [26], the depth-ratio values 
are consistent with the so-called ellipsoid shape (typical values 
comprised between 0.5 and 0.66), considered a common feature 
for volcanic lake basins, even though Rio Cuarto morphometry 
tends to approximate a steep-sided frustum model, corresponding 
to steep sides and flat bottom [27]. Such morphological features 
tend to prevent water vertical mixing, favoring meromictic 
conditions [28]. Thus, these physical parameters have a strong 
influence on the vertical distribution of chemical species, espeiially 
approaching the lake bottom where bio-geochemical processes 
have their maximum efficiency [29]. 

2.2 LImnological features and rollover events 

At Lake Hule, the limit between epi- and hypolimnion, marked 
by a very weak thermocline and the complete depletion of O2, was 
reported to occur at a depth ranging between —10 and — 12 m 
[23,24]. As reported by [22], this lake shows a persistent vertical 
stratification and the presence of COj in the deepest water strata. 
Occurrence offish death episodes, associated with sudden changes 
of water color from dark blue to red and strong smell in the lake 
surroundings, were reported by the local population in the last 
years (4 to 5 events from 1989 to 2002). These events, which took 
place during the cool, rainy and windy season (i.e. from December 
to February), were interpreted as caused by rollover phenomena 
[16,17,18,30]. 

The transition between epUimnion and hypolimnion in the 
meromictic Lake Rio Cuarto was measured at 20 and 25 m depth 
in May-June and January-February, respectively [18,25]. Rollover 
events, testified by fish killing and color changes of lake water from 
green to yellow-reddish, were observed in 1920 [31], between 
1978 and 1991 [22], in January 1997 [16] and in Fel)ruar)' 2010 
[18], just one month before our sampling. These events were 
possibly triggered by cooling of the shallow water layer caused by 
an anomalous weather characterized by low air temperature and 
strong winds [18,25,32]. 

Materials and Methods 

3.1 Sampling of water and dissolved gases 

Water and dissolved gas sampling was carried out in March 
2010 along vertical profiles from the lake surface to the bottom at 
regular intervals of 5 m (Lake Hule) and 10 m (Lake Rio Cuarto), 
in sites corresponding to the deepest points. Permission to sample 
in both lakes was guaranteed by Red Sismologica Nacional and 
Universidad de Costa Rica. According to the single hose method 
[33-35], water and dissolved gas samples were collected using a 
sampling line consisting of 10 m long Rilsan tubes (<1> = 6 mm) 
connected among them by steel connectors. Once the tube end 
was lowered to the chosen depth, water was pumped up to the 
surface through the sampling line using a 150 mL glass syringe 
equipped with a three-way teflon valve and transferred into plastic 
bottles after the displacement of a water volume double than the 
inner volume of the tube. One filtered (0.45 (tm) and two filtered- 
acidified (with ultrapure HCl and HNO3, respectively) water 
samples were collected in polyethylene botdes for the analysis of 
anions, cations and trace species, respectively. A fourth water 
ahquot was collected in glass bottles with the addition of HgCl2 for 
the analysis of water isotopes and '"^C/'^C ratios of total dissolved 
inorganic carbon (TDIC). Five hundred mL of water were filtered 
immediately after the sampling recovery through sterile cellulose 
mixed esters 0.22 |im pore size filters (GSWP, MiUipore, USA) for 
the analysis of prokaryotic populations. The filters were stored at 
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Figure 1. Map of Costa Rica witli the location of IHule and Rio Cuarto lalces. Modified after Alvarado et al. [18]. 

doi:1 0.1 371/journal.pone.01 02456.g001 



— 20°C in RNAlater solution (Quiagen, Italy), to prevent nucleic 
acid degradation. Dissolved gases were sampled using pre- 
evacuated 250 mL glass vials equipped with a Teflon stopcock 
and connected to the sampling line used to collect water samples. 
Sampling flasks were filled with water up to Vt of the inner volume 
[36-38]. 

3.2 Field measurements 

Water depth (m), temperature (°C), pH, Eh and electrical 
conductivity (EC; |J,S cm ') along the lake vertical profiles were 
measured using a Hydrolab MiniSonde 5 equipped with a data 
logger for data storage. The nominal precisions were: depth 
±0.05 m; T±0.rC; pH±0.2; Eh±20 mV; EC±1 nS cm"\ 
Alkalinity was measured in situ by acidimetric titration using 0.01 
N HCl. The analytical error for alkalinity analysis was S5%. 



3.3 Chemical and isotopic analysis of water and dissolved 
gases 

Main anions (CI , S04^ , NO3 , Br and F ) and cations 
(Na"", K"", Ca^+, Mg^+, NH/ and Li+) were analyzed by ion- 
chromatography (IC) using Metrohm 761 and Metrohm 861 
chromatographs, respectively. The analytical error for major water 
constituents was S5%. Trace elements at selected depths were 
analyzed at the INGV of Palermo by Inductively Coupled Plasma 
Mass spectrometry (ICP-MS, Agilent 7500-ce). For most of the 
elements the analytical uncertainty was in the order of 5-10% [39]. 

The "'0/"'0 and ^H/ H isotopic ratios of water (expressed as 
S'^O-HsO and 5D-H2O %o vs. V-SMOW, respectively) from 
selected depths were analyzed using a Finnigan Delta plusXP 
continuous-flow mass spei:trometer (MS) coupled with a Gas- 
benchll gas-chromatographic device (GBII), according to equili- 
bration techniques with CO2 for oxygen [40], and with for 
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Figure 2. Panoramic view and bathymetric map of Lal<e Hule (modified after Gdcice [24]). 

doi:1 0.1 371/journal.pone.01 02456.g002 
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Figure 4. Vertical profiles of temperature (°C, a), electrical conductivity (EC, in \iS cm~\ b), pH (c), and redox potential (Eh, in mV, d) 
in Lake Hule (blue line) and Lake Rio Cuarto (red line). 
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hydrogen [4-1]. The analytical uncertainties were ±0.08%o and 
±l%o for S'^O and 5D values, respectively. 

The '^C'''2C ratios of TDIC (expressed as 5''^Ctdic %o vs. V- 
PDB) at selected depths were determined on CO2 produced by 
reaction of 3 mL of water with 2 mL of anhydrous phosphoric 
acid in vacuum [42] using a Finningan Delta Plus XL mass 
spectrometer. The recovered CO2 was analyzed after a two-step 
extraction and purification procedures of the gas mixtures by using 
liquid N2 and a sohd-liquid mixture of liquid N2 and trichloro- 
ethylene [43,44]. The analytical uncertainty was ±0.05 %o. 

Dissolved gas composition was calculated using i) the compo- 
sition of the gas phase stored in the headspace of the sampling 
glass flasks, ii) the gas pressure in the flask headspace, iii) the 
headspace volume, and iv) the solubility coefficients in water of 
each gas compound [45] . The inorganic gas compounds hosted in 
the flask headspace (CO2, N2, CH4, Ar, O2, Ne, H2 and He) were 
determined using a gas-chromatograph (Shimadzu 15a) equipped 
with a Thermal Conductivity Detector (TCD). Methane was 
analyzed with a Shimadzu 1 4a gas-chromatograph equipped with 
a Flame Ionization Detector (FID). The analytical error for 
dissolved gas analysis was ^5%. 



The analysis of the '''C/'^C ratios of CO2 (expressed as S'^'C- 
CO2 %o vs. V-PDB) stored in the flask headspace (8'''C-C02strip) 
of selected samples was carried out with a Finningan Delta S mass 
spectrometer after standard extraction and purification procedures 
of the gas mixtures [43,44]. Internal (Carrara and San Vincenzo 
marbles) and international (NBS18 and NBS19) standards were 
used for the estimation of external precision. The analytical 
uncertainty was ±0.05%o. The ''^C/'^C ratio of dissolved CO2 
(S'-'C-COa) was calculated from the 5 ' C-CO2STRIP values using 
the Si factor for gas-water isotope equilibrium proposed by Zhang 
et al. [46], as follows: 



CO2 



C02s™p = (0.0049 X r)- 1.31 (1) 



The analysis of the ''^C/'^C and ^H/'H ratios of dissolved CH4 
(expressed as 5''^C-CH4 %o vs. V-PDB and 5D-CH4 %o vs. V- 
SMOW, respectively) of selected samples was carried out by mass 
spectrometry (Varian MAT 250) according to the procedure and 
the sample preparation described by SchoeU [47]. The analytical 
uncertainty was ±0.1 5%o. 
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The '^He/*He ratios, expressed as R/Ra values, where R is tlie 
^He/*He isotopic ratio in gas samples and Ra is that of the air 
equal to 1.39x10"'' [48,49], were determined in selected aas 
samples stored in the sampling flask headspace at the INGV 
laboratories of Palermo, using the method described in Inguag- 
giato and Rizzo [50]. The R/Ra values were corrected for air 
contamination on the basis of measured He/Ne ratios. The 
analytical uncertainty was ±0.3%. 

3.4 Microbiological analysis 

DNA extraction for the analysis of microbial populations was 
performed according to the protocol reported by MapeUi el al. 
[51] and quantified by NanoDrop 1000 spectrophotometer 
(Thermo Scientific, Waltham, MA). 16S rRNA gene was amplified 
in PGR reactions using universal primers for bacteria with GC- 
clamp as described in Marasco et al. [52]. Denaturing Gradient 
Gel Electrophoresis (DGGE), applied to the bacterial 16S rRNA 
gene amplified from the total water metagenome, was performed 
by loading DGGE-PGR products (-150 ng) in a 0.5 mm 
polyacrylamide gel (7% [w/v] acrylamide-bisacrylamide, 37.5:1) 
containing 40 to 55% urea-formamide denaturing gradient, where 
100% denaturant corresponds to 7 M urea and 40% [vol/vol] 
formamide [52]. DGGE profiles were analyzed by using Image J 
software (available at http://rsb.info.nih.gov/ij/) and cluster 
analysis was performed using Microsoft Excel XLSTAT software 
(Addinsoft Inc., New York, NY, USA). DGGE bands were excised 
from the gel, eluted in water, PGR amplified and sequenced as 
previously described [52]. The partial 16S rRNA gene sequences 
obtained from the excised DGGE bands were edited in Ghromas 
lite 2.01 (http://www.technelysium.com.au) and subjected to 
BLAST search (http://blast.ncbi.nlm.nih.gov/Blast.cgi). The nu- 
cleotide sequences were deposited in the EMBL pubhc database 
under the accession numbers HF930552-HF930593. To test the 
presence of bacteria involved in anaerobic ammonium oxidation 
(anammox), the functional gene hisA was amplified using primers 
hzsA_526F and hzsA_1857R as previously reported [53]. 



454 pyrosequencing assays were performed by using universal- 
bacterial primers targeting the variable regions of the 1 6S rRNA, 
V1-V3 (27 F mod 5' - AGRGTTTGATCMTGGCTCAG - 3'; 
519 R mod bio 5' - GTNTTACNGCGGCKGCTG - 3'), 
amplifying a fragment of approximately 400 bp, and 16S rRNA 
archaeal primers arch344F (5' - ACGGGGYGCAGCAGGCG- 
CGA- 3') and arch915R (5' - GTGCTCCCCCGCCAATTCCT 
-3'). The amplified 16S rRNA regions contained enough 
nucleotide variability to be useful in identification of bacterial 
and archaeal species [54,55] . PGR reactions and next generation 
454 pyrosequencing were performed at MR DNA laboratories 
(Shallowater, TX - U.S.A.). A first quality filtering was applied, 
removing all the sequences that were shorter than 300 bp, longer 
than 500 bp or with an average quality score under 25. All original 
and non-chimeric 454 sequences are archived at EBI European 
Read Archive. The high-quality 16S rRNA gene sequences 
obtained by 454 pyrosequencing were analysed using QIIME 
[56]. The sequences were clustered into operational taxonomic 
units based on a threshold of 97% (OTU97) sequence identity, 
using uclust [57] and one sequence for each OTU97, as 
representative, was aligned to Greengenes (http:/ /greengenes.lbl. 
gov/) using PyNast [56]. Sequence identification was conducted 
using Ribosomal Database Project classifier [58], with default 
parameters. For each sample rarefaction curves of the observed 
species and of Shannon index were estimated in order to analyse 
the species sampling coverage. The OTU97 diversity within and 
between sample/s (respectively alpha and beta diversity) was 
estimated using QIME workflow script alpha_rarefaction.py. 
Shannon diversity index was calculated by PAST software [59]. 
Library coverage was calculated for each library using the 
equation C= [1- (nl/N)] xlOO, where nl is the number of 
singleton OTU97, and N is the total number of reads in the library. 
To remove noise from the data, including potential rare 
contaminants, OTU97 not meeting the criterion of being present 
at least 0. 1 % of the total number of reads were removed. 
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Results 

4.1 Vertical profiles of temperature, EC, pH and Eh 

Temperature, EC, pH, and Eh along the vertical profiles of the 
lakes are shown in Tab. 1 and Fig. 4. Both Hule and Rio Cuarto 
lakes showed relatively high temperature at the surface (24. 1 and 
27.9°C, respectively), and a thermocline at shallow depths (starting 
from —2.5 and —5 m, respectively), with minimum temperatures 
of 20.8 and 24.6°C, respectively, at the lake bottoms (Fig. 4a). The 
temperature profiles were consistent with those reported in 
previous studies [17,18,23,24,25,32,60], except those of the 
epUimnion, likely because present and past measurements were 
carried out in different periods of the year. Lake Hule did not show 
a clear chemocline, as shown by the EC values that almost 
constantiy increased (from 84 to 140 |aS cm ) with depth (Fig. 4b). 
Conversely, Lake Rio Cuarto showed two chemoclines: the first 
one (from 159 to 186 |J,S cm ) near the surface and the second 
one (from 190 to 378 \xS cm"') between —40 and — 67 m depth. 
The vertical profile of pH values at Lake Hule exhibited a sharp 
decrease from 6.9 to 6.3 between the depths of 0 m and 10 m, and 
an opposite trend below this depth, where pH rose from 6.3 to 6.6 
(Fig. 4c). At Lake Rio Cuarto the pH values decreased in the 
shallower water strata (from 7.5 to 6.8) and from —40 to —60 m 
depth (from 6.8 to 6.5), and slightly increased (up to 6.6) at the lake 
bottom (Fig. 4c). Eh values at Lake Hule (Fig. 4d) showed a sharp 
decrease between —10 and — 15 m (from 33 to —163 mV) and 
reached the minimum values at lake bottom (—200 mV), whereas 
at Lake Rio Cuarto it strongly decreased (fi"om +166 at surface to 
— 191 mV) at the depth of 10 m displaying the lowest value (— 
246 mV) at the lake bottom. 

4.2 Chemical and isotopic composition of water samples 

Both lakes showed low TDS values (up to 90 and 246 mg L 
respectively, at lakes bottom) and a Ca^"'"-HC03 composition 
(Tab. 1). Concentrations of HCO,", NH4"^, Fe,„, and Mn (Fig. 5a- 
b) tended to increase towards the two lakes bottom (up to 61 and 
179 mg L"', 0.3 and 11 mg L"', 8 and 22 mg L"', 0.9 and 
0.7 mg L ' in Hule and Rio Cuarto, respectively), whilst oxidized 
nutrients NOa" and S04^~, typical electron acceptors in 
anaerobic environments, showed an opposite behaviour in Lake 
Rio Cuarto, decreasing to 0.03 and 0.4 mg L , respectively 
(Fig. 5b). On the contrary, F", CI", Ca^"", Mg-'"*", Na"", K"" and, 
only in Lake Hule, NO3 and S04^ , did not display specific 
vertical trends along the lakes water column. 

The SD-HjO values in Hule and Rio Cuarto lakes ranged from 
-26.5 to -20.5 %o and -24.4 to -19.7 %o V-SMOW, 
respectively, while those of 5"'0-H20 varied from —4.7 to — 
4.6 %o and from -4.5 to -4.0 %o V-SMOW, respectively (Tab. 1). 
The S'^'CxDic values were between —14.3 and —11.8 %o and — 
8.6 to —3.7 %o "V-PDB, in Hule and Rio Cuarto, respectively. 

Trace element composition did not differ significantly between 
the two lakes. The most abundant trace elements (>4 |j,g L ') 
along Hule and Rio Cuarto vertical profiles were Al, B, Ba, Rb, Sr 
and Zn. The maximum concentrations of Co, Cu, Ni, Ti and V (< 
2.2 |J.g L" ') were observed at the bottom layer of Lake Rio Cuarto 
(—67 m) and the other measured trace elements (As, Cd, Cr, Cs, 
Li, Mo, Pb, Sb, Se, Th, U) were aU <1 |ig L"' (Tab. 2). In terms 
of vertical distribution, those trace elements that clearly increased 
towards both lakes bottom were Al, As, Ba, Co, Ni, Sr, Ti and V 
(Tab. 2), whilst Mo concentrations showed a decrease with depth 
only in Lake Rio Cuarto. 
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Figure 6. Vertical profiles (in ^mol L M of CO2, N2, CH4 and O2 

doi:1 0.1 371 /journal.pone.OI 02456.g006 

4.3 Chemical and isotopic composition of dissolved gases 

Molecular nitrogen was the most abundant dissolved gas in the 
shallow portion of the two lakes (down to the depths of — 1 5 m and 
— 20 m at Lake Hule and Lake Rio Cuarto, respectively; Tab. 3). 
At lower depths CO, dominated the gas composition (up to 1090 
and 2090 |tmol L~' at Lake Hule and Lake Rio Cuarto, 
respectively), except at the bottom of Lake Rio Cuarto (Fig. 6a- 
b) where CH4 concentrations up to 2830 |J,mol L ' were 
measured. O2 is not present below — 10 m depth at Hule and 
Rio Cuarto, defining a clear anaerobic zone (Fig. 6a-b). Ar and 
Ne did not vary significantly with depth, whereas H2 and He 
concentrations increased with depth in both lakes (up to 0.01 and 
0.03 |imol L~' and to 0.04 and 0.3 |.tmol L~' in Hule and Rio 
Cuarto, respectively; Tab. 3). It is noteworthy to point out that He 
was an order of magnitude more abundant at Rio Cuarto than at 
Hule. The maximum total pressure (pTOT; Tab. 3) value of 
dissolved gases was measured at the bottom of Lake Rio Cuarto 
(2.9 atm), whereas pTOT in Lake Hule ranged from 0.79 to 1.1 
atm. 

The 5 ■ C-CO2 value at the bottom of Lake Hule was — 16.2 %o 
V-PDB (Tab. 3). At Lake Rio Cuarto, the S'-^C-COj values 
showed an increase with depth, ranging from — 14.3 at —20 m to 
—6.5 %o V-PDB at the lake bottom. No specific trends were 
recognized in the epHimnion (Fig. 7). The 5''*C-CH4 values, 
basically characterized by the same interval (from —77.2 to —60.7 
%o V-PDB) in both lakes, showed a rapid decrease in the Rio 
Cuarto hypolimnion. The 5D-CH4 values of Lake Rio Cuarto 
were significantly more negative (from —251 to —233 %o V- 
SMOW) when compared to that of Lake Hule bottom (—159 %o 
V-SMOW; Tab. 3). The R/Ra values, corrected for the presence 
of atmospheric helium [61], were 0.95 in Lake Hule (lake bottom) 
and 1.15 and 1.09 in Lake Rio Cuarto (at —50 and —67 m depth, 
respectively; Tab. 3). 




|jmol L 

in Lal<e Hule (a) and Lake Rfo Cuarto (b). 



4.4 Prokaryotic diversity along the water column 

Phylogenetic analyses of 16S rRNA DOGE derived sequences 
(Fig. 8a-b) allowed to detect 7 phyla within the bacterial 
communities and to identify the prevalent taxonomic groups 
colonizing the Hule and Rio Cuarto lakes at different depths 
(Tab. 4). Overall, the sequences were related to uncultured 
unclassified bacteria previously described in aquatic environments, 
mainly represented by freshwater of lacustrine origin. 

At Lake Hule a clear shift in taxa distribution was evaluated, 
corresponding to the transition at ~10 m depth of the redox 
potential from positive to negative. The lake epUimnion was 
mainly colonized by aerobic heterotrophic Bacteroidetes and 
Betaproteobacteric while deeper anoxic layers (>10 m depth; 
Fig. 4d) were inhabited by bacteria belonging to the phylum 
Chlorobi, comprising anaerobic photoautotrophic bacteria [Chlor- 
ohium clathratiforme and Ignavibacterium album). Bacteroidetes 
and Betaproteobacteria phyla were also the main components of 
the bacterial community in Lake Rio Cuarto. In this lake the 
shallower portion (down to the depth of 40 m) was colonized by 
Cyanobacteria afiiliated to the genera Synechococcus, Merismope- 
dia and Cyanobium. Differently from Lake Hule, the more 
uniform composition of the bacterial community in Lake Rio 
Cuarto can be related to the homogeneity of the redox conditions 
along the water column, which is negative in all the analyzed 
layers except at the lake surface (Fig. 4d). 

The results of DGGE analysis were taken into account to select 
a sub-set of samples to gain a deeper insight into the microbiome 
structure by massive pyrosequencing of bacterial and archaeal 16S 
rRNA libraries. This high-throughput analysis was applied to 3 
samples for each lake (0, 10, 15m depth from Lake Hule, named 
HO, HIO and HI 5, and 30, 50, 60 m depth from Lake Rio 
Cuarto, named RC30, RC50 and RC60). Unfortunately, any 
archaeal library could not be obtained from sample HO. The 
number of final reads varied among the samples, similarly to the 
OTU97 number, nonetheless a significant coverage of bacterial 
and archaeal diversity was reached in all the samples (Tab. 5). The 
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Figure 7. Vertical distribution of 8"C-C02 and S'^C-CH, of Lake Rio Cuarto. See the text for further details. 
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number of OTUgy present in the archaeal communities was 
constant along the water column of Lake Rio Cuarto, while in 
Lake Hule a significant increase was observed with depth (Tab. 5). 
In all the samples, Proteobacteria were the most abundant 
bacterial phylum, with the exception of the water samples 
collected from Lake Rio Cuarto at 50 and 60 m depths (RC50 
and RC60) where Cyanobacteria and Chloroflexi were the 
prevalent phyla, respectively (Tab. 6). Cyanobacteria were also 
present at high percentage (29.4%) in the oxic surface water 
sample in lake Hule (Tab. 6). The phylum Chlorobi was 
widespread in both the lakes in all the samples characterized by 
negative Eh values, with significant prevalence at 10 and 15 m 
depth in Lake Hule (18.5 and 17.6%, respectively). Among 
Proteobacteria, the EpsUon-subgroup was a minor community 
component in both lakes and Deltaproteobacteria were more 
abundant in Rio Cuarto, especially in the deeper layers (Tab. 6). 
Alpha- and Gamma-proteobacteria were differendy distributed in 
the two lakes. The latter were particularly abundant in shallower 
Hule layers (HID and HI 5), while the former were present at high 
percentages throughout the whole Hule water column (Tab. 6). 
The class Betaproteobacteria, mainly represented by the Coma- 
monadaceae and Methylophilaceae families, was abundant at all 
depths in both the lakes (Tab. 6). In Lake Hule between 12.9 and 



22.8% of the bacterial community was represented by sequences 
belonging to the ACK-Ml cluster of the order Actinomycetales, 
whose presence in lacustrine habitats was previously reported 
(Tab. 6) [62]. At the oxic-anoxic interfaces, anaerobic ammonium 
oxidation (anammox) was indicated as an autotrophic denitrifica- 
tion metabolism co-responsible of nitrogen loss from water 
environments [63] . The research of bacterial taxa known to be 
responsible of anammox reaction was performed by amplifying 
with specific primers the functional gene hzsA, encoding for 
hydrazine synthase and recently proposed as an anammox 
phylomarker [53]. The PCR amplification showed negative 
results, confirming that anammox populations are absent at Hule 
and Rio Cuarto lakes. 

As far as the archaeal community is concerned, Euryarchaeota 
were the most abundant phylum in Lake Rio Cuarto (up to 99%). 
Methanomicrobia were the most abundant class within this 
phylum, encompassing in particular the orders Methanomicro- 
biales and Methanosarcinales (Tab. 7). Lake Hule showed a 
different archaeal community, being dominated by Parvarchea 
and Micrarchaea, with significant concentrations of Crenarch- 
aeota (8.1 and 13.7% at 10 and 15 m depth, respectively), and a 
minor percentage of Methanomicrobia and unknown taxa (Tab. 7). 
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Figure 8. DGGE analysis performed on tlie bacterial 16S rRNA gene, showing the structure of the bacterial community inhabiting 
freshwater samples collected from the Hule and Rio Cuarto lakes (a); taxonomic identification of bacterial 16S rRNA sequences 
excised from DGGE bands cut from the Lake Hule and Rio Cuarto water profiles (b). 

doi:l 0.1 371/journal.pone.01 02456.g008 



Discussion 

5.1 Processes controlling the water chemistry 

Water isotopes can provide notable information on physical- 
chemical processes afTecting the chemistry of volcanic lakes, such 
as evaporation, water-rock interaction and hydrothermal/ meteor- 
ic inputs [64]. As shown in Fig. 9, water samples plot nc'ar the 
Global Meteoric ^Vater Line (GMWL) [6.5] and the Costa Rica 
Surface Water Line [66], indicated that in both lakes the main 
water source is meteoric, consistently with their Ca^"'"(Mg^"*)- 
HCO3 composition, which is typical for superficial waters and 
shallow aquifers worldwide [67]. Both lakes show a slight D- and 
"'O- depletion at increasing depth, likely related to evaporation 
affecting epUimnetic waters [64,68,69]. 

The parallel increases of HCO3 (Fig. 5) and dissolved CO2 
(Fig. 6) along the vertical profiles suggest that the behaviour of 
these two chemical species is controlled by the following reactions: 

C02+H20^H2COi (2) 

and 



H2COi+H20^HiO+ + HCOi " (3) 

The observed weak decreases of S04^~ and NO^" concentra- 
tions (Fig. 5) with depth possibly result from microbial activity 
occurring at anaerobic conditions. The lack of free oxygen in the 
hypoUmnion is favorable for nitrate reduction by microbial 
denitrification, a typical process in anoxic water bodies [70-74]. 
The genus Pseudomonas is known to include denitrifier species 
[7.'>] and was retrieved at high abundance in pyrosequencing 
libraries in the anoxic layers of Lake Hule, constituting up to 16'% 
of the total bacterial community (Tab. 6). In the Hule anoxic 
layers, 16S rRNA pyrosequencing allowed to detect additional 
denitrifying genera like Sulfuricurvum, Opitutus and Geothrix 
(Tab. 6). Sulfate reducing bacteria (SRB) of the genus Syntropho- 
bacter were retrieved by 16S rRNA pyrosequencing in the deepest 
layers of the Rio Cuarto water column (Tab. 6), and could be 
responsible of the weak depletion observed for S04^ (Fig. 5b). 
Nevertheless, the relatively low S04^ and NO3 concentrations, 
typical of meteoric-sourced lakes, implies that sulfate reduction 
and denitrification have a minor impact on the chemistry of the 
two lakes. The increase of NH4 concentrations with depth (Fig. 5) 
is apparendy suggesting direct NH4^ production within the 
hypoUmnion via ammonification processes [76]. 

The increase of Fe and Mn contents in the deepest water layers 
can be attributed to direct production inside the bottom sediments 
by minerogenic processes [77-79], although their presence as 
solutes is limited by the formation of insoluble Fe- and Mn- 
hydroxides. Gocke [24] suggested that the high concentration of 
Fe in Lake Hule is also caused by the addiction of yellow/ 
brownish Fe(OH)3 material through the southern brooklets, which 
subsequently precipitates in the hypolimnion and iron is reduced 
to the ferrous state, as also supported by the relatively low Eh 
values (Fig. 4). Oxidation of hypolimnetic Fe^""" in the epilimnion 
would explain the yellow-reddish color of the shallow water layer 
that was occasionally observed in these lakes as a consequence of 



water rollover [18,25]. Nevertheless, the red coloration observed 
at Lake Hule in February 1991 was likely caused by the presence 
of dense purple clumps or masses floating of Merismopedia [18], a 
genus belonging to the phylum Cyanobacteria that were observed 
by DGGE at -30 and -40 m depth in Lake Rio Cuarto (Tab. 4). 

As shown in the spider-diagrams of Fig. 1 0, where concentra- 
tions of Al, Ba, Cr, Cu, Ni, Rb, Sr, Ti and V at maximum depths 
for both lakes are normalized to those measured in basalt rock 
samples collected from the young intra-caldera cone at Laguna 
Hule (the only one available) [80], water-rock interactions 
efficiently mobilized soluble elements such as Ba, Rb and Sr, 
whereas Al and Ti were basically retained in the rock matrix. In 
particular, Cr and Ni, as well as As and Co, are possibly related to 
the dissolution of Mn-and Fe-oxide particles that setded through 
the chemocline [78,81,82]. The concentrations of dissolved V are 
strongly correlated with those of Fe, similarly to what observed for 
Mo and Mn [29,83], likely because they belong to the same 
mineralogical paragenesis. For what concerns the other trace 
elements, Cu and Zn may be related to dissolution of stable 
organic complexes buried in the bottom sediments [29]. Cs, Rb 
and B, which are strongly correlated with Li (Tab. 2), can be 
considered as conservative elements, likely due to the strong 
affinity of alkali ions and boric acid for the aqueous phase [82]. 
The relatively low Mo concentrations at increasing depth in Lake 
Rio Cuarto (Tab. 2) may be related to its consumption during 
microbial nitrate reduction [29]. 

5.2 Processes governing chemical and isotopic 
composition of dissolved gases 

5.2.1 Noble gases, N2, O2, and H2. Dissolved gas species in 
volcanic lakes basically originate from i) external sources (e.g. 
atmosphere, volcanic-hydrothermal fluids) and/or ii) microbial 
activity occurring both in lake water and at water-sediment 

interface [4,14,28,84,85,86]. 

Dissolved Ar and Ne in lakes are related to air dissolution 
through the lake surface, a process that is mainly controlled by 
atmospheric pressure and the water temperature [87]. The inert 
noble gases behave inerdy in any bio-geochemical process and 
thus along the lake water column they are affected by advection 
and diffusion. Accordingly, Ar and Ne concentrations in the two 
investigated lakes did not show significant variations with depth 
(Tab. 3). Conversely, O2, which is typically consumed by aerobic 
microbial populations for oxidation of organic matter and reduced 
ionic species, rapidly decreases with depth, and was virtually 
absent at depths >5 and 10 m, in Hule and Rio Cuarto lakes, 
respectively. It is worth noting that the N2/Ar ratios were slighdy 
higher than that of air saturated water (~40), suggesting the 
addition of N2 from an extra-atmospheric source. This hypothesis 
is expected to be confirmed by 6'^N values that are presentiy not 
available, although the relatively high N2/Ar ratios are apparendy 
consistent with nitrate depletion with depth and microbial 
denitrification in both lakes. Consistendy with the N2 excess, the 
distribution of N2 concentrations in both lakes showed significant 
variations with depth (Tab. 3), probably related to N2 production 
and consumption by denitrifiers and nitrogen fixing prokaryotes, 
respectively. Microbial N2 fixation, depending on light [88] and 
the presence of bio-avaUable trace metals [89], can be carried out 
by heterocyst-forming species in water and in sediment pores [90- 
92]. Cyanobacteria were indeed retrieved by both DGGE and 
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pyrosequeiicing in surface layers of Rio Cuarto and Hule lakes 
(where they constitute 26% of the total bacterial community in the 
oxic layer HO, Tab. 6), supporting the occurrence of N2 frxation in 
both the lakes. 

H2 increase with depth in the hypolimnion at Hule and Rio 
Cuarto (Tab. 3) suggests a production of H2 likely related to 
fermentation of organic matter under anaerobic conditions at the 
water-sediment interface. Additionally, photoreactions carried out 
by Cyanobacteria, abundantly present in the Rio Cuarto deep 
layers and in the upper layer of the Hule water columns (Tab. 6), 
could be responsible of Hj production [93-97]. Once produced at 
the lake bottom, H2 can be consumed acting as electron donor for 
hydrogeno trophic methanogenic archaea and SRB [98-100], 
detected in Rio Cuarto pyrosequencing libraries. Moreover, it 
slowly diffuses up to shallower, oxygenated layers where it can be 
consumed by hydrogen-oxidizing bacteria [101-104]. 

The presence of an extra-atmospheric source for helium can be 
recognized on the basis of the R/Ra values (Tab. 3), which are 
relatively high (up to 20 or more) for mantle gases, and as low as 
0.01 in fluids from crustal sources [61]. Dissolved gas samples from 
Hule and Rio Cuarto lakes showed R/Ra values ~ 1 that, coupled 
with the relatively high He/Ne ratios (49 and 4.1 at Lake Rio 
Cuarto and Lake Hule, respectively), indicate a significant fraction 
of mantie He, whose uprising is likely favored by the fault system 
characterizing this area [18]. 

5.2.2 CO2 and CH4. CO2 and CH4 are the most abundant 
extra-atmospheric dissolved gases present in Hule and Rio Cuarto 
lakes. As already mentioned, dissolved CO2 controls pH values 
and HC03~ concentrations. Previous studies [17,18,20,22,23] 
have hypothesized that these lakes are alfected by CO2 inputs 
through the bottom, as supported by the presence of C02-rich 
bubbling pools and caverns or boreholes with high CO2 
concentrations characterizing this area [18,105]. A significant 
contribution of mantle CO2 is indicated by the 5'''C-C02 value of 
the dissolved gas sample collected at the maximum depth of Lake 
Rio Cuarto (—6.6 %o vs. V-PDB; Tab. 3), which is in the range of 
mantie gases (from -8 to -4 %o vs. V-PDB) [106]. Although not 
confirmed by the 5'"'*C-C02 values, the CO2/CH4 ratio measured 
in the dissolved gas at the bottom of Lake Hule (4.7) is too high, 
even higher than that of Rio Cuarto bottom sample (0.63), to be 
ascribable to microbiological processes. This would imply that 
even at Lake Hule a strongly negative isotopic signature of CO2 is 
externally added to the bottom waters, possibly from a C02-rich 
source deriving from oxidation of previously produced hydrocar- 
bons. 

The 5''^C-C02 values at the bottom of Lake Hule (—16.2 %o vs. 
V-PDB) and at depths between —20 and —50 m in Lake Rio 
Cuarto (as low as — 14.3 %o vs. V-PDB; Tab. 3) were intermediate 
between those generated by organic matter degradation [24] and 
mantie degassing [107-109], indicating that along the vertical 
profiles of both lakes, excluding the bottom layers, biogenic 
processes are the most important sources of CO2. 

According to the classification proposed by Whiticar [110], the 
5''^C-CH4 and 5D-CH4 values of the Hule and Rio Cuarto lakes 
indicate that CH4 has a biogenic origin (Fig. 1 1). The vertical 
profiles of the concentrations and 5' '^C values of CO2 and CH4 of 
Lake Rio Cuarto (Fig. 7) were thus produced by the combination 
of different processes occurring at various depths in the lake: 

1) At the bottom of the lake, CO2 inputs from a deep source 
likely related to the hydrothermal fluid circulation [18,111] 
promote methanogenic processes that have their maximum 
efficiency within the sediments. Methanogenesis takes place 
through i) CO2 reduction and ii) degradation of organic 
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Table 5. Library coverage estimations and sequence diversity of 16S rRNA. 





Sample 


N. reads/sample 


N. OTU,7 


% Coverage* 


Shannon index** 


HO Bacteria 


9384 


586 


0.98 


4.21 


H10 Bacteria 


11115 


615 


0.98 


4.21 


HIS Bacteria 


1S872 


1260 


0.97 


4.98 


RC30 Bacteria 


32932 


3017 


0.95 


4.60 


RC50 Bacteria 


13291 


1609 


0.94 


4.37 


RC60 Bacteria 


13S30 


1882 


0.93 


5.08 


HIO Arclnaea 


140S 


68 


0.99 


2.76 


HIS Arcliaea 


5889 


289 


0.98 


2.08 


RC30 Archaea 


2429 


177 


0.96 


2.66 


RCSO Archaea 


2005 


172 


0.95 


2.81 


RC60 Archaea 


2937 


178 


0.97 


2.34 



^Library coverage was calculated as C = 1-n/N, where n is the number of OTU97 without a replicate, and N is the total number of sequences. 
**Shannon diversity Index calculated using PAST. 
doi:1 0.1 371 /journal.pone.Ol 02456.t005 



matter through acetate fermentation [47,110,112-115]. 
These processes can be described by the following reactions: 

C02 + ?>H+ +?,e- ^CH^ + H20 (4) 

and 

* CH^COOH^ * CH^ + CO2 (5) 

where the * indicates the intact transfer of the methyl position 
to CH4. 

2) In the hypolimnion, microbial CH4 production is stiU active, 
although this process is accompanied by CO2 dissolution, 
CH4 oxidation, and vertical diffusion of both the gas species. 
Moreover, in correspondence of aerobic/ anaerobic boundar- 
ies, anaerobic decomposition of organic matter [116-118], 
and CH4 oxidation carried out by methanotrophic bacteria 
can efficiently produce CO2 in lakes [86,1 19-121]. 

3) In the epilimnion, photosynthetic microorganisms (e.g. 
Cyanobacteria) convert light into biochemical energy through 
oxygenic photochemical reactions combined with CO2 
assimilative reduction. Vertical water circulation favors the 
activity of photosynthetic and methanotrophic bacterial 
populations, as well as the continuous addition of atmospheric 
gases from the lake surface.. 

These hypotheses were confirmed by the 16S rRNA pyrose- 
quencing of samples collected along the water column of Lake Rio 
Cuarto, demonstrating that archaeal communities encompass 
almost exclusively methanogenic populations (Tab. 7) typical of 
freshwater ecosystems, namely Methanomicrohiales and Methano- 
sarcinales [122-124], as also observed in freshwater meromictic 
lake sediments [125]. Methanosarcinales included solely the 
acetate-utilizing methanogen Methanosaeta, the most abundant 
archaeal genus along the Rio Cuarto water column. Within the 
H2-CO2 utilizing methanogens of the order Methanomicrohiales, 
Methanoregula was the prevalent genus, but unclassified Metha- 
nomicrohiales and Methanoregulaceae sequences were also detect- 
ed (Tab. 7). 

The lack of isotopic data along the vertical profile of Lake Hule 
did not allow to investigate in detail the (bio)-geochemical 



processes controlling the vertical profiles of CO2 and CH4. In 
this lake the majority of the archaeal 16S rRNA sequences were 
affiliated within unclassified Euryarchaeota, showing high similar- 
ity with the Candidate divisions Micrarchaea and Parvarchaea 
(Tab. 7) previously described by metagenomics studies of an acidic 
ecosystem by Baker et al. [126,127]. These archaeal sequences 
belong to the ARMAN (Archaeal Richmond Mine Acidophilic 
Nanoorganisms) lineages, which are among the smallest cellular 
life forms known [126], still poorly described from an ecological 
perspective. The presence of novel uncultivated lineages in the 
Lake Hule water is linked to neither specific metabolism nor the 
influence on the water and dissolved gas chemistry. However, 
besides a minor fraction of known acetotrophic methanogenic 
Methanosarcinales (Tab. 7), the archaeal community of Lake Hule 
included also the Miscellaneous Crenarchaeota Group (MCG), 
within the phylum Crenarchaeota (Tab. 7). MCG is a cosmopol- 
itan clade that was previously detected in both freshwater [128] 
and marine ecosystems [129], where it had been hypothesized to 
have a significant role in dissimHatory methane oxidation [129]. 
This hypothesis leads to the speculation that MCG could have the 
same ecological function also in the Lake Hule. It is worth noting 
that the minor percentage of known methanogenic archaea in 
Lake Hule compared to that of Lake Rio Cuarto corresponds to 
the differences between the lakes in CH4 concentrations (Tab. 3). 

16S rRNA pyrosequencing of bacterial communities showed 
that type I and type II methanotrophic bacteria, belonging to the 
Gamma- (i.e. Methylocaldum, Methylomonas, Crenothrix) and 
Alpha-subgroup of proteobacteria (i.e. Methylocystaceae) 
[125,130], respectively, were abundant in the anoxic layers of 
Hule and Rio Cuarto (Tab. 6), suggesting a key role in the carbon 
cycle. Within the Beta-proteobacteria, additional families that 
encompass methylotrophic bacteria, namely Methylophilaceae , 
Rhodocyclaceae , and Comamonadaceae [131,132], were retrieved 
by deep sequencing in the same water layers both in Lake Hule 
and Lake Rio Cuarto, the latter hosting up to 36% of 
Methylophilaceae at 30 m depth (Tab. 6). Within the family 
Comamonadaceae, relevant in Lake Hule, 5.2% of the bacterial 
sequences from the surface layer were affiliated to the genus 
Limnohahitans , which was reported to play a functional key role in 
freshwater habitats and showing high ecological diversification 
[133]. Moreover, 6.3% of the bacterial sequences were affiliated to 
the genus Ruhrivivax that includes, among the few characterized 



PLOS ONE I www.plosone.org 



17 



July 2014 I Volume 9 | Issue 7 | e102456 



Bio-Activity Volcanic Lakes: Hule and Rio Cuarto (Costa Rica) 





o 

LO 


o 
o 


o 
o 




0^ 


o 


LO 


0^ 


o 
o 


rN 

00 




o 
o 






m 

CO 


LO 
U-) 


9.52 


1^ 
o 




in 


o 
rN 




CO 


c^ 
m 


0^ 
O 


o 


0^ 
O 


>0 
ro 


o 
o 


o 
o 


O 
O 


o 
o 


o 
o 


!^ 


o 


o 


d 


d 


d 


d 






d 


d 


d 


d 


d 


d 


u-i 




ro 


d 




d 




d 


d 


nn 


d 


d 


d 


d 


d 


d 


d 


d 


d 


a\ 
o 


m 


o 
o 


o 
o 


CO 


o 


Lfl 

Ln 


o 


00 


o 
o 


o 


in 


o 
o 


CO 

o 


rM 


LO 


00 

ro 


O 


a\ 

00 


o 


o 


o 
o 


Ln 




m 
p 




o 




ro 




rM 
O 


O 
O 


o 
o 


m 
o 


\6 


o 


d 


d 


d 


n-i 


d 


d 




d 


d 




d 


d 


d 


rn 


d 








d 


d 


d 


d 




d 


d 


d 


d 


d 


d 


d 


d 


d 


00 




o 


o 


in 




fN 






o 


Pi 




o 




ro 


rv! 






p 


CTv 




o 


ro 




o 


? 


o 




So 




o 


o 


o 


o 


rn 


o 


d 


d 


rn 


fN 


d 


d 


ro 


d 


d 




d 


d 


d 






■=3-" 






d 


d 


d 


d 




d 


d 


d 


d 


d 


d 


d 


d 


d 




CO 


in 


CO 


m 


q 


o 




O 


o 


o 


O 




CT; 
O 


m 




00 


O 
p 






o 
o 


o 
o 


o 










o 
^' 


o 






m 




fN 


o 


o 




rn 


d 




d 


d 


d 


d 


d 


d 


d 






rN 






d 




d 


d 


d 


d 


d 


d 






d 


d 


d 


d 


d 




o 
o 




ro 




CO 

o 


O 

fN 


o 
o 


o 
o 


o 


o 
o 


o 
o 


O 
o 




O 


rn 


rv. 


o 


O 


o 
o 


in 


o 
o 


o 
o 


o 
o 


o 


m 




1^ 


LO 

0^ 


o 
o 


CO 

o 


CM 


LO 






o 


o 


d 




d 




d 


d 


d 


d 


d 


d 


d 






rN 


d 


d 


d 




d 


d 


d 


d 


d 


d 


d 




d 


d 


d 


d 


d 




^ 


o 
o 


o 
o 


o 
o 


u-i 


rv. 


o 
o 


o 
o 


o 
o 


ro 


ro 
O 


o 
o 


rN 


o 
o 


o 
o 


o 
o 


o 
o 


o 
o 


o 


oi 


o 
o 


o 
o 


o 
o 




LO 


o 


o 


ro 

CO 


o 
o 


ro 




CO 


LO 


m 
o 


o 


o 


d 


d 


d 




d 


d 


d 


ro 


d 


d 




d 


d 


d 


d 


d 


d 




d 


d 


d 


d 


d 


d 


d 




d 


d 


ro 


Lfi 




d 



E S 



I 1 o I g 



ri t; 



cu u tl u 
o y 



1-1 Ti "D P 



u u u c ^ ^ 



E 
E 



"a _• _: 



re fp n 
-tr: o 



E E 



tl tl 



^ ^ -d o 

O O =) X 



-D 

■ — ^-j +-> 



u u ^ u u 



0 _IJ _u _y _u _u 



cQ cQ CD ca 



S o 



~ ~ ~ ~ 



^-' S2 



E £ E 



On; H; lU 

i-i 4-1 i-i ^-l ^-l ^-l 



4-' -M O 



CO CO CO 



PLOS ONE I www.plosone.org 



18 



July 2014 I Volume 9 | Issue 7 | e102456 



Bio-Activity Volcanic Lakes: Hule and Rio Cuarto (Costa Rica) 



o 








yo 


o 


o 


o 


U 


o 


o 


o 


cc 


d 


d 


d 



o 
u 

n 







































o 




o 






o 




o 


o 


o 


ro 


o 


o 


o 


o 


o 






d 


d 


d 


d 




d 


d 


d 


d 


d 




d 


d 


d 


d 


d 


d 


d 


m 






LO 






IN 


o 




o 




o 










fN 


fN 






o 
















1^ 




o 


o 












d 


d 


rn 


d 






d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 






1^ 


u-i 




u-i 






Lfl 


o 




o 


o 




o 


o 










p 




0^ 










p 




p 


p 


p 


p 


p 


Ln 












































































vO 




CO 




O 


vO 


O 


LO 














o 




ro 


O 


rn 


q 


a\ 


p 


p 


CO 


p 


p 


p 


ro 








p 


in 


p 


p 


p 
























































LO 


















LO 


o 


o 




O 


o 


o 




LO 




'=r 




m 








o 


o 


p 


o 


o 


o 


o 


o 


o 




o 


i< 






VD 


m 


p 


CO 


o 


o 




d 


d 


d 


d 


d 


d 


d 


d 






d 




d 




d 


d 


d 


O 


o 


o 


o 


o 


o 


o 




o 


o 


1^ 






m 


Oi 




o 


o 


O 


o 


o 




o 


o 


o 




o 


o 


o 


o 


o 


CO 


CO 




o 


o 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 




d 


d 



o 



5 o 



u O u 



!3= y u ti 



^ -C 



U t-J - I— 



;2 u u 0 



O "d "u ^ ~ ^ 



5 E 



JS o 



2 b 



-ii ^ J:; o 



E E -Si 



^ .i 



™ o 



"t^ ^ 



en cn CD ca 



qj tu OJ OJ o y 
-t; -t; >^ n 









fD 


fD 


fD 


fD 


fD 


OJ 


OJ 


OJ 




OJ 


OJ 


a 


OJ 


u 


U 


U 


U 


u 


U 


u 


u 


fD 


TO 


fD 


TO 


fD 


fD 


fD 


fD 




X! 


X! 








JD 




o 


O 


O 


O 


o 


O 


O 


o 


OJ 


OJ 


OJ 




tu 


OJ 


a 


OJ 


o 


O 


o 


o 


O 


o 


o 


o 


Q. 


D. 


Q. 


Q. 


Q. 


Q. 


Q. 


Q. 


fD 


rs 


fD 


nj 


fD 


fD 


fD 


fD 


E 


E 


E 


E 


E 


E 


E 


E 


E 


E 


E 


E 


E 


E 


E 


E 








nj 


fD 


fD 


fD 


fD 








e? 



























n I "JJ ^ "JJ 
UJ +j 4_< +J 

t-^ U U U 



O +j +j 



cr _: 

. fD 

"a £ 



E E £ 



s_ j_ ro 



PLOS ONE I www.plosone.org 



19 



July 2014 I Volume 9 | Issue 7 | e102456 



Bio-Activity Volcanic Lakes: Hule and Rio Cuarto (Costa Rica) 



u 



m 
I 



o 
I 



3 
Z 



^ fij (ij 



vo CO a\ 



o .y u < 
^ ^ 5 > 



species, strains able to oxidize carbon monoxide producing carbon 
dioxide and hydrogen [134]. The presence of the genus 
Syntrophobacter at 60 m depth in Rio Cuarto (RG60) is in 
agreement with the establishment in deep anoxic layers of 
syntrophic relations between organic acid degrading bacteria 
and methanogenic archaea. Members of this genus were 
commonly detected in anaerobic mixed cultures, where they 
obtain energy from the anaerobic oxidation of acetate, growing 
syntrophically with hydrogen- and formate-utilizing methanogenic 
archaea [135]. The RC60 sample showed a high percentage of 
sequences affiliated to the order Dehalococcoidetes (Tab. 6), which 
comprises obligate organohalide respirers, widely detected in 
marine and freshwater ecosystems [136,137]. The presence of 
organohalide compounds favors the competition with methano- 
gens for the use of molecular hydrogen [138]. Hence the finding of 
Dehalococcoidetes in the deeper layers of Lake Rio Cuarto, 
retrieved by both pyrosequencing (RC50) and DGGE (RC60), 
suggest the presence of naturally occurring organo-halogens in the 
water that could serve as electron acceptors for organohalide- 
respiring bacteria. 

Further confirmation of the importance of anaerobic microbial 
processes on the CO2-CH4 balance can be obtained by comparing 
measured S'^^Ctdic values with those expected assuming isotopic 
equilibrium between CO2 and HCO3 . Isotopic fractionation 
caused by the reaction between dissolved CO2 and HCOj" is 
quantified by the enrichment factor (S2), as follows [139]: 



S2 = S^^C-HC0i- -S^^C-C02 = 94S2/T(K)-23,.9 (6) 



by: 



Theoretical S'^'Cxdic; values (S'^'CxDiCcaic) can be computed 



< 



^ fij n n 



-TDlCcalc- 



^13 C - CO2 + 82 X (7/CO3- ) / [(//CO3- ) + (CO2) 
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As shown in Fig. 12, water samples from the shallower strata 
(down to 40 m depth) of Lake Rio Cuarto displayed 5'^Ctdig and 
5''^CxDic:caic values basically consistent. On the contrary, samples 
from depth >40 m showed a strong difference between the two 
sets of values: at — 50 m depth, 5' '^CxDic<:ak: were more negative 
than 5''^CtdiC) whereas an opposite behavior was observed in the 
deeper water layer, as well as at the maximum depth of Lake Hule 
(Tab. 1). At the lake bottoms, continuous inputs of hydrothermal 
CO2, characterized by 5'''C-C02 values significantly less negative 
with respect to that already present in the lake, are likely 
responsible of the positive shift of the 5''^CxDiCcaic values, since 
this external CO2 was not in equihbrium with HCO3 . In the 
shallower layers, especially at the depth of —60 m, addition of 
non-equilibrated biogenic CO2 played an opposite role (Fig. 12), 
whereas at depth £40 m CO2 concentrations were too low to 
significantly affect the S'^^CxDiCcaic values, which were consistent 
with the S'^'CxDic ones. The disagreement between measured and 
calculated S'^'Cxdicj values, depending on both microbial activity 
and inputs of hydrothermal CO2, was documented in other 
meromictic lakes hosted in volcanic environments, such as Lake 
Kivu, D.R.C. [34] and the Italian lakes of Albano, Averno and 
Monticchio [86]. 

Although the multidisciplinary approach applied in the present 
study allowed to link the presence of different prokaryotic 
taxonomic groups to the observed physical conditions and the 
concentrations of chemical species along the water columns, the 
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Figure 9. 6^"0-5D diagram for the water samples from Lake Hule (blue squares) and Lake Rio Cuarto (red squares). See the text for details. 
doi:1 0.1 371/journal.pone.01 02456.g009 
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Figure 10. Spider-diagrams, where concentrations of selected trace elements in Lake Hule (a) and Lake Rio Cuarto (b) maximum 
depths are normalized to those measured in basalt rock samples collected from the young intra-caldera cone at Laguna Hule [80]. 
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Figure 11. S C-CH4 vs. SD-CH4 plot (modified after Whiticar [110]) of Lake Hule (blue square) and Lake Ri'o Cuarto (red squares). 

See the text for further details. 
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ecological function of certain prokaryotcs in these intriguing 
ecosystems, particularly in the case of Lake Hule, remains cryptic. 
In particular, among the bacterial community. Lake Hule hosted 
the Actinomycetales ACK-Ml cluster [140], whose phenotypic 
and metabolic traits have not yet been described. The ACK-Ml 
cluster was one of the most abundant bacterial taxonomic groups 
in Lake Hule, reaching up to 22.8% in the oxic water layer HO 
(Tab. 6). Moreover, in the Lake Hule waters, the Alphaproteo- 
bacterial order Rickettsiales showed relatively high concentrations 
(18.3% of the total bacterial community in the oxic layer HO; 
Tab. 6). This order comprises intracellular organisms, pinpointing 
the importance of symbiotic relationships in these lakes. In this 
context, the impact of the associations between bacteria and algae 
[141] or phytoplancton [142] on nutrients re-minerahzation was 
recently discussed showing the crucial role of trophic levels 



interaction on the food web of lacustrine habitats, possibly relevant 
also in volcanic lakes. 

Conclusions 

Hule and Rio Cuarto are meromictic maar lakes mainly fed by 
meteoric water, and characterized by significant amounts of 
dissolved gases, partially consisting of CO2 having a hydrothermal- 
magmatic origin, in their hypolimnion. They are currently 
classified as low activity or, alternatively, "Nyos-type" lakes [4], 
implying that a limnic eruption could be expected to occur from 
these lakes, as confirmed by the rollover events they have 
experienced. However, gases stored in the deep layers of Hule 
and Rio Cuarto are fundamentally different with respect to those 
of Nyos and Monovm lakes, a difference that must be considered 
for evaluating the eruption risk. The gas reservoirs of the two 
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Figure 12. Vertical distribution of measured and calculated S Ctdic for the water samples from Lake Rio Cuarto. See the text for 
further details. 
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Cameroonian killer lakes are composed of almost pure COj and 
basically their temporal evolution only depends on a high 
magmatic gas input rate [12,13]. At Nyos, the risk of gas bursts 
was successfully mitigated artificially by discharging the deep- 
seated gases at the lake surface [35,143]. On the contrary, the gas 
reservoirs of Hule and Rio Cuarto lakes consist of CO2, CH4 and 
N2 in comparable amounts, mainly controlled by the activity of a 
microbial network governed by CO2 and CH4 metabolism, thus 
the possible occurrence of a lake roUover that may pose a local risk 
is not directly related to the input rate of external CO2. 

Despite geographic separation. Lake Rio Cuarto and Lake Hule 
showed similar physical-chemical settings, though hosting phylo- 
genetically distinct bacterial and archaeal communities. Phyloge- 
netic difference apart, however, both lakes have revealed the 
presence of the same prokaryotic ecological functions deeply 
involved in affecting water and gas chemistry. 

On the whole. Lake Hule and Lake Rio Cuarto host a 
C02(CH4, N2)-rich gas reservoir which is mainly controlled by the 
complex and delicate interactions occurring between geosphere 
and biosphere and whose monitoring can appropriately be carried 
out by coupling the conventional geochemical approach with 



studies about prokaryotic colonization. Consequently, for these 
lakes we can introduce the new definition of hio-activity lakes. This 
term can be extended to several other volcanic lakes which show 
similar compositional features of water and dissolved gases, e.g. 
Kivu (D.R.C. -Rwanda) [34,144], Monticchio, Albano and Averno 
(Italy) [37,86,145-147], Pavin (France) [121,148]. 

Acknowledgments 

We wish to thank Lorenzo Brusca and Sergio BeUomo (INGV-Palermo) for 
their laboratory assistance for trace elements analyses. The authors would 
like to thank Corentin Caudron {Earth Observatory of Singapore) for the 
detailed and constructive reviews of the original manuscript. 

Author Contributions 

Conceived and designed the experiments: JC FT FM SB S. Calabrese DR 
GC RM BC RA OV GP S. Caliro RMA. Performed the experiments: JC 
FT FM SB S. Calabrese RM BC RA OV GP FC GB S. Caliro. Analyzed 
the data: JC FT FM SB. Contributed reagents/materials/ analysis tools: FT 
SB S. Calabrese DR GC OV GP CR RMA. Contributed to the writing of 
the manuscript: JC FT FM SB S. Calabrese DR RM BC OV GP RMA. 



PLOS ONE I www.plosone.org 



23 



July 2014 I Volume 9 | Issue 7 | e102456 



Bio-Activity Volcanic Lakes: Hule and Rio Cuarto {Costa Rica) 



References 

1. Rouwct D, Tassi F, Mora-Amador R, Sandri L, Chiarini V (2014) Past, present 
and future ot'volcanie lake monitoring. J Voleanol (jeotherm Res 272: 78-97. 

2. Brown G, Rymer H, DowdenJ, Kapadia P, Stevenson D, ct al. (1989) Energy 
budget analysis for Pods crater lake: implications for predicting volcanic 
activity. Nature 339: 370-373. 

3. Brantley SL, Agustsdottir AM, Rowe GL (1993) Grater lakes reveal volcanic 
heat and volatile fluxes. Geol Sec Am 3: 175-178. 

4. Pasternaek GB, Varekamp JC (1997) Volcanic lake systematics I. Physical 
constraints. Bull Voleanol 58(7): 528-538. doi:10.1007/s004450050160. 

5. Anzidei M, Carapezza ML, Esposito A, Giordano G, LeUi M, et al. (2008) The 
Albano Maar Lake high resolution bathymetry and dissolved C()2 budget 
(CoUi Albani volcano, Italy): constrains to hazard evaluation. J Voleanol 
Geotherm Res 171: 258-268. 

6. Hurst T, Christenson B, Gole-Baker J (2012) Use of a weather buoy to derive 
improved heat and mass balance parameters for Ruapehu Crater Lake. 
J Voleanol Geotherm Res 235: 23-28. 

7. Rouwet D, Tassi F (2011) Geochemical monitoring of volcanic lakes. A 
generalized box model for active crater lakes. Arm Geophys 54: 161—173. doi: 
10.4401/ag-5035. 

8. Klmg (;\V; Glark MA, Compton HR, De\-incJD, Evans VVC, ct al. (1987) The 
1986 Lake Nyos gas disaster in Cameroon, West Africa. Science 236: 169—175. 

9. Sigurdsson H, Devince JD, Tchoua FM, Presser TS, Pringle MKW, et al. 
(1987) (Origin of the lethal gas burst from Lake Monoun, Cameroon. J Voleanol 
Geotherm Res 31: 1-16. 

10. Barberi F, Ghelini W, MarineUi G, Martini M (1989) The gas cloud of Lake 
Nyos (Cameroon, 1986): Results of the Italian technical mission. J Voleanol 
Geodicrm Res 39: 125 131. 

1 1. Giggcnbach WF (1990) Water and gas chemistry of Lake Nyos and its bearing 
on the eruptive process. J Voleanol Cieotherm Res 42: 337-362. 

12. Evans WC, Kling GW, Tutde ML, Tanyileke G, White LD (1993) Gas buildup 
in Lake Nyos, Cameroon: the recharge process Euid its consequences. Appl 
Geochem 8: 207-221. 

13. Evans WC, White LD, Tutde ML, Kling GW, Tanyileke G, et al. (1994) Six 
years of changes at Lake Nyos, Cameroon, yield clues to the past and cautions 
for the future. Geochem J 28: 139-162. 

14. Kusakabc M (1996) flazardous crater lakes. In: Scarpa R, Tilling Rl, 
editors. Monitoring and mitigation of volcano hazards. Springe r-Vcrlag, Berlin, 
pp. 573-598. 

15. Rice A (2000) Rollover in volcanic crater lakes: a possible cause for Lake Nyos 
type disasters. J Voleanol Geotherm Res 97: 233-239. 

16. Haberyan KA, Horn SP, Umana GV (2003) Basic limnology of fifty-one lakes 
in Costa Rica. Rev Biol Trop 51: 107-122. 

17. Tassi F, VascUi O, Fernandez E, Duarte E, Martinez M, et al. (2009b) 
Morphological and geochemical features of crater lakes in Costa Rica: an 
overview. J Limnol 68: 193—205. 

18. Alvarado GE, Soto GJ, Salani FM, Ruiz P, Hurtado de Mendoza L (201 1) The 
formation ajid evolution of Hule and Rio CuEuto maars, Costa Rica. J Voleanol 
Geodierm Res 201: 342-356. 

19. Horn SP, Haberyan KA (1993) Physical and chemical properties of Costa 
Rican lakes. Natl Geogr Res Explor 9(1): 86-103. 

20. Horn SP (2001) The age of the Laguna Hule explosion crater, Costa Rica, and 
the timing of subsequent tephra eruptions: evidence from lake sediments. Rev 
Geol Am Cent 24: 57-66. 

21. Umana (i, Haberyan KA, Horn SP (1999) Limnology in Costa Rica. In: Gopal 
B, Wetzel RW, editors. Limnology in Developing Countries 2: 33—62. 

22. HaberyEui KA, Horn SP (1999) Chemical and physical characteristics of seven 
volcanic lakes in Costa Rica. Brenesia 51: 85-95. 

23. Umana G (1993) The planktonic community of Laguna Hule, Costa Rica. Rev 
Biol Trop 41(3): 499-507. 

24. Gocke K (1997) Basic morphometric and limnological properities of Laguna 
Hule, a caldera lake in Costa Rica. Rev Biol Trop 44/45: 537-548. 

25. Gocke K, Bussing W, Cortes J (1987) Morphometric and basic limnological 
properties of the Laguna de Rio Cuarto, Costa Rica. Rev Biol Trop 35(2): 277- 
285. 

26. Carpenter SR (1983) Lake geometry: implications for production and sediment 
accretion rates. J Theor Biol 105: 273-286. 

27. LehmanJT (1975) Reconstructing the rate of accumulation of lake sediment. 
The effect of sediment focusing. Quatem Res 5: 541-550. 

28. Martini M, Giannini L, Prari F, Tassi F, Capaccioni B, et al. (1994) Chemical 
characters of crater lakes in the Azores and Italy: the anomaly of the Lake 
Albano. GeochemJ 28: 173-184. 

29. Wetzel RG (2001) Limnology: Lake and River Ecosystems. 3rd Ed., Academic, 
San Diego, Calif., USA. 

30. Soto GJ (1999) Gcologia Regional de la hoja Pods (1: 50.000). In: Alvarado GE, 
Madrigal LA, editors. Estudio Geologico-Geot^cnico de Avance a factibihdad 
del P. Laguna Hule. Inf. Interne ICE, San Jose, Costa Rica. pp. 15-^5. 

31. Sapper K (1925) Los Volcanes de la America Central. Max Niemayer, Halle 
(Saale). 144 p. 

32. Gocke K, Bussing W, CortcsJ (1990) The annual cycle of primary productivity 
in Laguna de Rio Cuarto, a volcanic lake (maar) in Costa Rica. Rev Biol Trop 
38(2B): 387-394. 



33. Tassi F, Vaselli O, Giannini L, Tedesco D, Nencetti A, et al. (2004) A low-cost 
and effective method to collect water and gas samples from stratified crater 
lakes: the 485 m deep lake Kivu (DRC). Proc. lAVCEl Gen. Ass., Puchon, 
Chile, 14-19 November 2004. 

34. Tassi F, Vaselli O, Tedesco D, Montegrossi G, Darrah T, et al. (2009a) Water 
and gas chemistry at Lake Kivu (DRC): geochemicEil evidence of vertical and 
horizontal heterogeneities in a multi-basin structure Geochem. Geophys. 
Geosyst. 10, doi:10.1029/2008GC002191 

35. Tassi F, Rouwet D (2014) An overview of the structure, hazards, and methods 
of investigation of Nyos-type lakes from the geochemical perspective. J Limnol 
73(1): DOI: 10.4081/jlimnol.2014.836 

36. Chiodini G (1996) (iases dissolved in groundwaters: analytical methods and 
examples of applications in central Italy. In: Marini L, OttoncUo (}, editors. 
Proc. Symp. Environmental Geochemistry. Castelnuovo di Porto, Rome, 22—26 
May 1996. pp. 135-148. 

37. Caliro S, Chiodini G, Izzo G, Minopoh C, Signorini A, et al. (2008) 
Geochemical and biochemical evidence of lake overturn and fish kill at hake 
Avemo, Italy. J Voleanol Geotherm Res 178: 305-316. 

38. Tassi F, Vaselli O, Luchetti G, Montegrossi G, Minissale A (2008)Metodo per 
la determinazione dei gas disciolti in acque naturali. Int Rep CNR-IGG, 
Florence, n" 10450. 11 p. 

39. Calabrese S, Aiuppa A, Allard P, Bagnato E, Bellomo S, ct al. (2011) 
Atmospheric sources and sinks of volcanogenic elements in a basaltic volcano 
(Etna, Italy), (ieochim Cosmochim Acta 75: 7401—7425. 

40. Epstein S, Mayeda TK (1953) Variation of the O/ O ratio in natural waters. 
Geochim Cosmochim Acta 4: 213-224. 

41. Nelson ST (2000) A simple, practical metiiodology for routine VSMOW/SLAP 
normalization of water samples analyzed by continuous flow methods. Rapid 
Commun Mass Spectrom 14: 1044—1046. 

42. Salata GG, Roelke LA, Cifuentes LA (2000) A rapid and precise method for 
measuring stable carbon isotope ratios of dissolved inorganic carbon. Mar 
Chem 69: 153-161. 

43. Evans WC, White LD, Rapp JB (1998) (}eochcmistr\' of some gases in 
hydrothermal fluids from the southernjuan de Fuca ridge. J (ieophys Res 15: 
305-313. 

44. Vaselli O, Tassi F, Montegrossi G, Capaccioni B, Giannini L (2006) Sampling 
and analysis of fumju-oHc gases. Acta Vulcanol 1-2: 65-76. 

45. Whitfield M (1978) Activity coefiicients in natural waters. In: Pytkowicz RM, 
editor.Activity Coefficients in Electrolyte Solutions. CRC Press, Boca Raton, 
Florida, pp. 153-300. 

46. Zhang J, Quay PD, Wilbur DO (1995) Carbon isotope fractionation during 
gas-water exchange and dissolution of C02. Geochim Cosmochim Acta 59: 
107-114. 

47. Schoell M (1980) The hydrogen and carbon isotopie composition of methane 
from natural gases of various origins. Geochim Cosmochim Acta 44: 649—661. 

48. Mamyrin BA, Tolstikhin IN (1984) Helium isotopes in nature. Elsevier, 
Amsterdam. 

49. Ozima M, Podosek FA (2002) Noble Gas Geochemistry. Cambridge University 
Press, UK. 

50. Inguaggiato S, Rizzo A (2004) Dissolved helium isotope ratios in ground- 
waters: a new technique based on gas-water re-equilibration and its application 
to Stromboli volcanic system, -^pl Geochem 19: 665-673. http://dx.doi.org/ 
10.1016/j.apgeochem.2003.10.009 

51. Mapelli F, Varela MM, Barbato M, Alvarino R, Fusi M, et al. (2013) 
Biogeography of planktonic microbial communities across the whole Mediter- 
ranean Sea. Ocean Sci Discuss 10: 291-319. doi:10.5194/osd-10-291-2013 

52. Marasco R, RoUi E, Lttoumi B, Vigani G, Mapelli F, et al. (2012) A drought 
resistance-promoting microbiome is selected by root system under desert 
farming. PLoS ONE 7(10): e48479. doi:10.1371/journal.pone.0048479 

53. Harhangi HR, Le Roy M, van Alen T, Hu B-I, Groen J, et al. (2012) 
Hydrazine synthase, a unique phylomjuker with which to study the presence 
and biodiversity of anammox bacteria. Appl Environ Microbiol 78: 752-758. 

54. Van de Peer Y, Chapelle S, De Wachter R (1996) A quantitative map of 
nucleotide substitution rates in bacterial rRNA. Nucleic Acids Res 24(17): 
3381-3391. 

55. Chakravorty S, Helb D, Burday M, Connell N, Alland D (2007) A detailed 
analysis of 16S ribosomal RNA gene sc^gments for the diagnosis of pathogenic 
bacteria. J Microbiol Methods 69: 330- 339. 

56. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, et al. 
(2010) QIIME allows analysis of high-throughput community sequencing data. 
Nat Methods 7: 335-336. 

57. Edgar RC (2010) Search and clustering orders of magnitude faster than 
BLAST. Bioinformatics 26(19): 2460-2461. 

58. Wang Q, Garrity GM, Tiedje JM, ColeJR (2007) Naive Bayesian Classifier for 
Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy. 
Appl Environ Microbiol 73(16): 5261-5267. 

59. Hammer 0, Harper DAT, Ryan PD (2001) PAST: palcontological statistics 
sofhvare package for education and data analysis. Palaeontol Electronica 4(4): 
1-9. 



PLOS ONE I www.plosone.org 



24 



July 2014 I Volume 9 | Issue 7 | e102456 



Bio-Activity Volcanic Lakes: Hule and Rio Cuarto {Costa Rica) 



60. Umana Q (2010) Comparison of basic limnological aspects of some crater lakes 
m the Cordillera Volcanica Central, Costa Rica. Rev Geol Amer Central 43: 
137-145. 

61. Craig H, Lupton JE (1976) Primordial neon, helium and hydrogen in oceanic 
basalts. Earth Planet Sci Lett 31: 369-385. 

62. Wu QL, Zwart G, Schauer M, Kamst-van Agterveld MP, Hahn MW (2006) 
Baeterioplankton Community Composition along a Salinity Gradient of 
Sixlt-cn High-Mountain Lakes Located on the Tibetan Plateau, China. AEM 
72: 5478 5485. 

63. Zhu G, Jetten MSM, Kusehk P, Ettwig KF, Yin C (2010) Potential roles of 
anaerobic ammonium and methane oxidation in the nitrogen cycle of wetland 
ecosystems. Appl Microbiol Biotechnol 86: 1043-1055. 

64. Varekamp JC, Kreulen R (2000) The stable isotope geochemistry of volcanic 
lakes, with examples from Indonesia. J Volcanol Geotherm Res 97: 309-327. 

65. Craig H (1961) Isotopic variations in meteoric waters. Science 133: 1702-1703. 

66. Laehniet MS, Patterson WP (2002) Stable isotope values of Costa Rican surface 
waters. J Hydrol 260: 135-150. 

67. Berner EK, Berner Rj\ (1987) Global Water Cycle: (ieochemistrv and 
Environment. Prentice-Hall, Inc, Englewood Clilfs, Newjersey. p. 397. 

68. Matsubaya O, Sakai H (1978) D/H and ^^'O/^^O fractionation factors in 
evaporation of water at 60 and 80°C. GeochemJ 12: 121-126. 

69. Rowe GLJr (1994) Oxygen, hydrogen and sulfur isotope systematics of the 
crater lake system of Poas volcano, Costa Rica. GeochemJ 28: 263—287. 

70. Alexander M (1961) Introduction to Soil Microbiology. John Wiley & Sons, 
New York. p. 472. 

71. Buresh RJ, Patrick WH (1981) Nitrate reduction to ammonium and organic 
nitrogen in an estuarine sediment. Soil Biol Biochem 13: 279-283. 

72. Stewart WDP, Preston T, Peterson HG, Christofi N (1982) Nitrogen cycling in 
cutrophic frcshwaters. Philosoph Transact Royal Soc B 296: 491—509. 

73. Ahlgren I, Sorensson F, Waara T, Vrede K (1994) Nitrogen budgets in relation 
to microbial transformations in lakes. Ambio 23(6): 367—377. 

74. Brune A, Frenzel P, Cypionka H (2000) Life at the oxic-anoxic interface: 
microbial activities and adaptation. FEMS Microbiol Rev 24(5): 691-710. 

75. Carlson CA, Ingraham JL (1983) Compairison of denitrification by Pseudomo- 
nas stutzeri, Pseudomonas aeruginosa, and Pamcoccus denitrificans . Appl 
Environ Microbiol 45: 1247-1253. 

76. Molongoski JJ, Klug MJ (1980) Anaerobic metabolism of particulate organic 
matter in the sediments of a hypereutrophic lake. Freshwater Biol 10: 507-518. 

77. Davison W, Heaney SI, TallingJF, Rigg E (1980) Seasonal transformations and 
movements of iron in a productive English lake with deep water anoxia. 
Schweiz Z Hydrol 42: 196-224. 

78. Balistrieri LS, Murray JW, Paul B (1992) The cycling of iron and manganese in 
the water column of Lake Sammamish, Washington. Limnol Oceanogr 37: 
510-528. 

79. Hong\'e D (1997) Cycling of iron, manganese, and phosphate in a meromictic 
lake. Limnol Oceanogr 42: 635—647. 

80. Prosser JT, Carr MJ (1987) Poas volcano, Costa Rica: geology of the summit 
region and spatial and temporal variations among the most recent lavas. 
J Volcanol Geotherm Res 33: 131—146. 

81. Balistrieri LS, Murray JW, Paul B (1994) The geochemical cycling of trace 
elements in a biogenic meromictic lake. Geochim Cosmochim Acta 58(19): 
3993-4008. 

82. VioUier E, Jezequel D, Michard G, Pepe M, Sarazm G, et al. (1995) 
Geochemical study of a crater lake (Pavin Lake, Prance): trace-element 
behaviour in the monimolimnion. Chem (ieol 125;1 2): 61 72. 

83. Schallcr T, Moor HC, Wehrli B (1997) Reconstructing the iron cycle Irom the 
horizontal distribution of metals in the sediment of Baldeggersee. Aquat Sci 59: 
326-344. 

84. Varekamp JC, Pastemack GB, Rowe GLJr (2000) Volcanic lake systematics II. 
Chemical constraints. J Volcanol Geotherm Res 97: 161-179. 

85. Schmid M, Halbwachs M, Wehrli B, Wuest A (2005) Weak mixing in Lake 
Kivu: new insights indicate increasing risk of uncontrolled gas eruption. 
Geochem Geophys Geosyst 6: 1—11. 

86. Gabassi J, Tassi F, Vaselh O, Fiebig J, Nocentim M, ct al. (2013) 
Biogeochemical processes involving dissolved CO2 and CH4 at Albano, 
Averno, and Monticchio meromictic volcanic lakes (Central-Southern Italy). 
Bull Volcanol 75(1): 1-19. 

87. Weiss R (1970) The solubility of nitrogen, oxygen and argon in water and 
seawater. Deep Sea Res 17: 721-735. 

88. Tison DL, Palmer FE, StaleyJT (1977) Nitrogen fixation in lakes of the Lake 
Washington drainage basin. Water Res 11: 843-847. 

89. Hyenstrand P, Blomqvist P, Pettersson A (1998) Factors determining 
cyanobacterial success in aquatic systems - a literature review. Arch Hydrobiol 

15: 41-62. 

90. Moeller RE, Roskoski JP (1978) Nitrogen-frxation in the littoral benthos of an 
oligotrophic lake. Hydrobiologia 60(1): 13-16. 

91. Loeb SL, Renter JE (1981) The epilithic periphyton community: a five-lake 
comparative study of community productivity, nitrogen metabolism and depth- 
distribution of standing crop. Verh Intemat Verein Limnol 21: 346-352. 

92. Valiela I (1991) Ecology of coastal ecosystems. In: Barnes RSK, Mann KH, 
editors. Fundamentals of aquatic ecology. Blaekwell Science, Oxford, pp. 57— 
76. 

93. Benemann JR, Weare NM (1974) Hydrogen evolution by nitrogen-fixing 
Anabaena cylindrical cultures. Science 184: 174—175. 



94. (ireenbaum E (1982) Photosynthetic hydrogen and oxygen production: kinetic 
studies. Science 215: 291-293. 

95. Asada Y, Kawamura S (1986) Aerobic hydrogen accumulation by a nitrogen- 
fixing Cyanobacterium, Anabaena sp. Appl Environ Microbiol 51: 1063-1066. 

96. Asada Y, Miyake J (1999) Photobiological hydrogen production. J Biosci 
Bioengineer 88(1): 1-6. 

97. Bandyopadhyay B, StockelJ, Min H, Sherman LA, Pakrasi HB (2010) High 
rates of photobiological H2 production by a cyanobacterium under aerobic 
conditions. Nature Communications 1: 139. doi:10.1038/ncommsl 139 

98. Mah RA, Ward DM, Baresi L, Glass TL (1977) Biogenesis of methane. Annu 
Rev Microbiol 31: 309-341. 

99. Zehnder AJB (1978) Ecology of methane formation. In: Michell R, editor. 
Water pollution microbiologyk. J. WUey & Sons Inc, New York. pp. 349-376. 

100. Thauer RK, Badziong W (1980) Respiration with sulfate as electron acceptor. 
In: Knowles CJ, editor.Diversity of bacterial respiratory systems. CRC Press, 
Boca Raton, Fla, 2. pp. 65-85. 

101. Aragno M, Sehlegel HG (1981) The hydrogen-oxidizing bacteria. In: Starr MP, 
Stolp H, Friiper HG, Ballows A, Sehlegel HG, editors.The prokaryotes. A 
handbook of habitats, isolation and identification of bacteria. Vol. 1. Springer- 
Verlag, Berlin. 

102. Bowien B, Sehlegel HG (1981) Physiology and biochemistry of aerobic 
hydrogen-oxidizing bacteria. Ann Rev Microbiol 35: 405^52. 

103. Conrad R, Aragno M, Seiler W (1983) Production and consumption of 
hydrogen in a eutrophic lake. Appl Environ Microbiol 45: 502-510. 

104. Bianchi L, Mannelli F, Viti C, Adessi A, De Phifippis R (2010) Hydrogen- 
producing purple non-sulfur bacteria isolated from the trophic lake Avemo. Int 
J Hydr En 35: 12213-12223. 

105. Zimmer MM, Fisher 'FP, Hilton UH, Alvarado GE, Sharp ZD, et al. (2004) 
Nitrogen systematics and gas fluxes of subduction zones: insights from Costa 
Rica arc volatiles. Geochem Geophys Geosyst 5(5): 1-19. doi:10.1029/ 
2003GC000651 

106. Barnes I, Irwin WP, White DE (1978) Global distribution of carbon dioxide 
discharges and major zones of seismicity. US Geological Survey, Water- 
Resources Investigation, 78-39, Open File Report. 

107. O'Leary MH (1988) Carbon isotopes in photosynthesis. BioScience 38: 328- 
336. 

108. Rollinson H (1993) Using geochemical data: evaluation, presentation, 
interpretation. Longman Scientific and Technical, New York, p. 352. 

109. Hoefs J (2009) Stable Isotope Geochemistry, 6th edn. Springer, Berlin, 
Germany, p. 288. 

110. Whiticar MJ (1999) Carbon and hydrogen isotope systematics of bacterial 
formation and oxidation of methane. Chem Geol 161: 291—314. 

111. Alvarado GE, Soto GJ, Pullinger CR, Escobar R, Bonis S, et al. (2007) 
Volcanic activity, hazards, and monitoring. In: Bundscliuh J, Alvarado GE, 
editors. Central America: Geology, Resources and Hazards, Vol. 2. Taylor & 
Francis, London, pp. 1155-1188. 

1 12. Mah RA, Ward DM, Baresi L, Glass TL (1977) Biogenesis of methane. Annu 
Rev Microbiol 31: 309-341. 

113. Barker J F, Fritz P (1981) Carbon isotope fractionation during microbial 
mediane oxidation. Nature 293: 289-291. 

114. Schoell M (1988) Multiple origins of methane in the Earth. Chem Geol 71: 1- 
10. 

115. Woese CR, Kandler O, Wheelis ML (1990) Towards a natural s\stem of 
organisms: proposal for the domains Archaca, Bacteria and luiearva. Proe Natl 
Acad Sci 117: 1 1576 1 1579. 

1 16. RuddJWM, Hamilton RD, Campbell NLR (1974) Measurement of microbial 
oxidation of methane in lake water. Limnol Oceanogr 19: 519—524. 

117. Rich PH (1975) Benthic metabolism of a soft-water lake. Verh Internat Verein 
Limnol 19: 1023-1028. 

118. Rich PH (1980) Hypolimnetic metabolism in three Cape Cod lakes. Amer 
Midland Natur 104: 102-109. 

119. Frenzel P, Thebrath B, Conrad R (1990) Oxidation of methane in the oxic 
surface layer of a deep lake sediment (Lake Constance). FEMS Microbiol Eeol 
73: 149-158. 

120. Casper P (1992) Methane production in lakes of different trophic state. Arch 
Hydrobiol Beih Ergebn Limnol 37: 149-154. 

121. Lopes F, Viollier E, Thiam A, Michard G, Abril G, et ah (2011) 
Biogeochemical modeling of anaerobic vs. aerobic methane oxidation in a 
meromictic crater lake (Lake Pavin, France). Appl Geochem 26: 1919—1932. 

122. Franzmann PD, Liu YT, Balkwill DL, Aldrich HC, deMacario EC, et al. (1997) 
Methanogenium frigidum sp. nov., a psychrophilic, H2-using methanogen from 
Ace Lake, Antarctica. Int J Syst Bacteriol 47: 1068-1072. 

123. Brauer SL, Cadillo-Quiroz H, Ward RJ, Ya\itt JB, Zinder SH (2011) 
Methanoregula boonei gen. nov., sp. nov., an acidiphilie methanogen isolated 
from an acidic peat bog. IntJ Syst Evol Microbiol 61: 45-52. 

1 24. Chaudhary PP, Brableova L, Buriankova I, RuHk M (20 1 3) Molecular diversity 
and tools for deciphering the methanogen community structure and diversity in 
freshwater sediments. Appl Microbiol Biotechnol 97: 7553-7562. 

125. Borrel G, J6z6quel D, Biderre-Petit C, Morel-Desrosiers N, Morel J, et al. 
(201 1) Production and consumption of methane in freshwater lake ecosystems. 
Res Microbiol 162: 832-847. 

126. Baker BJ, Tyson (iW, Webb Rl, Flanagan J, Hugenholtz P, et al. (2006) 
Lineages of acidophilic Archaea revealed by community genomic analysis. 
Science 314: 1933-1935. doi:10.1 126/science.l 132690. 



PLCS ONE I www.plosone.org 



25 



July 2014 I Volume 9 | Issue 7 | e102456 



Bio-Activity Volcanic Lal<es: Hule and Rio Cuarto (Costa Rica) 



127. Baker BJ, ComoUi LR, Dicka GJ, Hauscr LJ, Hyatt D, ct al. (2010) Enigmatic, 
ultrasmall, uncultivated Archaea. PNAS 107: 8806-8811. 

128. Borrel G, Lehours A-C, Grouzet O, Jiziquel D, Rockne K, et al. (2012) 
Stratification of Archaea in the deep sediments of a freshwater meromictic lake: 
Vertical Shift fi-om Methanogenic to Uncultured Archaeal Lineages. PLoS 
ONE 7:e43346. doi:10.1371/journal.pone.0043346. 

129. Biddle JF, Lipp JS, Lever MA, Lloyd KG, Sorensen KB, et al. (2006) 
Heterotrophic Archaea dominate sedimentary subsurface ecosystems off Peru. 
PNAS 103: 3846-38,') 1. 

130. Stoeeker K, Bendinger B, Sehoning B, Nielsen PH, Nielsen JL, et al. (2006) 
Cohn's Crenothrix is a filamentous methane oxidizer with an unusual methane 
monooxygenase. PNAS 103: 2363-2367. 

131. Ghistoserdova L, Kalyuzhnaya MG, Lidstrom ME (2009) The expanding 
world of methylotrophie metabolism. Annu Rev Microbiol 63: 477-499. 
doi: 1 0. 1 1 46/annurev.micro.09 1 208.073600 

132. Beck DAC; Kalyuzhnaya MG. Malfatti S, IVinge SG, del Rio TO, et al. (2013) 
.\ melageiiomie insight into iresliwarer metliane-utilizing communities and 
evidence for cooperation between the Methylococcacme, and the Methylophi- 
laceae. PeerJ l:e23. doi:10.7717/pecrj.23 

133. JezberaJ, Jezberova J, Kasalicky V, Simek K, Hahn MW (2013) Patterns of 
Limnohabitans microdiversity across a large set of freshwater habitats as 
revealed by reverse line blot hybridization. PLoS ONE 8: e58527. doi:10.1371/ 
journal.pone. 0058527. 

134. Maness P, Huang J, Smolinski S, Tek V, Vanzin G (2005) Energy generation 
from the CO oxidation-hydrogen production pathway in Ruhrivivax 
gekltinosus. Appl Environ Microbiol 71: 2870-2874. 

135. Harmsen HJM, Van Kuijk BLM, Pluggc CM, Akkcrmans ADL, De Vos WM, 
et al. (1998) Syntrophobact^r fummoxidans sp nov, a syntrophic propionate- 
degrading sulfate-reducing bacterium. Int J Syst Bact 48: 1383-1387. 

136. Hug LA, GastcUe GJ, Wrighton KG, Thomas BG, Sharon I, et al. (2013) 
Community genomic analyses constrain the distribution of metabolic traits 
across the Chloroflexi phylum and indicate roles in sediment carbon cycling. 
Microbiomc 1:22, doi: 10.1 186/2049-2618-1-22. 

137. Zanaroli G, Balloi A, Negroni A, Borruso L, Dalfonchio D, et al. (2012) A 
Chloroflexi bacterium dechlorinates polychlorinated biphenyls in marine 



sediments under in situ-like biogeochemical conditions. J Haz Mat 209—210: 

449-457. 

138. Balloi A, Rolli E, Marasco R, Mapelli F, Tamagnini I, et al. (2010) The role of 
microorganisms in bioremediation and phytoremediation of polluted and 
stiessed soils. Agrochimica 54(6): 353-369. 

139. Mook WG, Bommerson JC, Staverman WH (1974) Carbon isotope 
fractionation between dissolved bicarbonate and gaseous carbon dioxide. 
Earth Planet Sci Lett 22: 169-176. 

140. Zwart G, Crump BC, Agterveld tkl, Hagen F, Han SK (2002) Typical 
freshwater bacteria: an aiialvsis of available 16S rRNA gene sequences from 
plankton of lakes and rivers. Aquat Microb Ecol 28: 141-155. 

141. Eigemann F, Hilt S, Salka I, Grossart H (2013) Bacterial community 
composition associated with freshwater algae: species specificity vs. dependency 
on environmental conditions and source community. FEMS Microbiol Ecol 83: 
650-663. 

142. Paver SF, Hayek KR, Gano KA, Fagcn JR, Brown CT, et al. (2013) 
Interactions between specific phytoplankton and bacteria affect lake bacterial 
community succession. Environ Microbiol 15: 2489-2504. 

143. Kusakabc M, Ohba T, Issa YY, Satake H, Ohizumi T, et al. (2008) Evolution 
of CO2 in lakes Monoun and Nyos, Cameroon, before and during controlled 
degassing. GeochemJ 42: 93—118. 

144. Schoell M, Tietze K, Schoberth SM (1988) Origin of methane in Lake Kivu 
(East-Central Africa). Chem Geol 71: 257-265. 

145. Carapezza ML, LcUi M, Tarchini L (2008) Geochemistry of the ,\lbano and 
Ncmi crater lakes in the volcanic district of Alban Hills (Rome, Italv). J Voleanol 
Geotherm Res 178: 297-304. 

146. Caracausi A, Nuccio PM, Favara R, Nicolosi M, Paternoster M (2009) Gas 
hazard assessment at the Monticchio crater lakes of Mt Vulture, a volcano in 
Soutiiern Italy. Terra Nova 21: 83-87. 

147. Chiodini G, Tassi F, Galiro S, Chiarabba C, Vaselli O, et al. (2012) Time- 
dependent CO2 variations in Lake Albano associated with seismic activity. Bull 
Voleanol 74: 861-871. 

148. Aeschbaeh-Hcrtig W, Hofer M, Kipfer R, Imboden DM, Wider R (1999) 
Accumulation of mande gases in a permanentiy stratified volcanic lake (Lac 
Pavin, France). G«ochim Cosmochim Acta 63: 3357-3372. 



PLOS ONE I www.plosone.org 



26 



July 2014 I Volume 9 | Issue 7 | e102456 



